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Abstract 

We present a study of the many-body interaction between a D* and multi-p. We use an extrapolation 
to SU(4) of the hidden gauge formalism, which produced dynamically the resonances /2(1270) in the pp 
interaction and 1)2(2460) in the pD* interaction. Then let a third particle, p, D* , or a resonance collide 
with them, evaluating the scattering amplitudes in terms of the Fixed Center Approximation of the Faddeev 
equations. We find several clear resonant structures above 2800 MeV in the multibody scattering amplitudes. 
They would correspond to new charmed resonances, D3, D\, and Z3g, which are not yet listed in the 
PDG, which would be analogous to the P3(1690), /4(2050), P5(2350), /6(2510) and /sr^lTSO), Kl{2M^), 
Kl{2'i%Q) described before as multi-p and i^*-multi-/9 states respectively. 



PACS numbers: 



I. INTRODUCTION 



One of the important aims in the study of the strong interaction is to understand the nature and 
structure of hadronic resonances. The search for new resonances is another goal both in theories 
and experiments. At low energy, using the input of chiral Lagrangians [l|-[5(] and implementing 
unitarity in coupled channels, one develops a theoretical tool, chiral unitarity theory, which explains 
the two-body interaction very successfully @-17|. For the three-body interaction, the pioneer work 
of Ref. 18] combined Faddeev equations and chiral dynamics and reported several S-wave = 
resonances which qualify as two mesons-one baryon molecular states. One more step, using the 
Fixed Center Approximation Ul^] to Faddeev equations for multi-/9(770) states, was given in 
Ref. fii] in which the resonances /2(1270), p3(1690)(3-~), /4(2050)(4++), /95(2350)(5~ ), and 
/6(2510)(6''~''') were explained as basically molecules of an increasing number of /3(770) particles 
with parallel spins. Similarly, it was also found in Ref. that the resonances -Rr|(1430), ^"1(1780), 
K|(2045), (2380) and a new could be understood as molecules made of an increasing number 

of p(770) and one K*(^92) meson. The A5 + (2000) puzzle was explained in Ref. 
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system with the same method 

The Fixed Center Approximation to Faddeev equations is technically very simple and fairly 
accurate when dealing with bound states, as discussed in Ref. [27|]. It should also be avoided when 
dealing with states which have enough energy to excite its components in intermediate states [30(]. 
It should be noted that the results of the Fixed Centre Approximation to Faddeev equations have 
been confirmed by the variational method approach with the effective one-channel Hamiltonian in 
a more recent work 3l[ which predicts a narrow DNN quasi-bound state in the energy range of 
about 3500 MeV. 

In our present work, the main ideas of Refs. 



23|,|2j] are followed to search i^*-multi-ff resonances 



in the charm sector. Usmg effective Lagrangia;;rof the local hidden gauge theory BH, the pp 
interaction was studied in Ref. 36|| with on-shell factorized Bethe-Salpeter equations. It was found 
that the pp interaction was attractive in the isospin zero, spin and 2 channels, particularly in the 
tensor channel where it led to the formation of a pp quasibound state or molecule that could be 
associated to the /2(1270) ( I{J^^) = 0(2++) ). With the same formalism, the composite systems 
of light {p and w) and heavy {D*) vector mesons were studied in Ref. [s^]. In that work, a strong 
attraction was found in the isospin, spin channels (/,£') = (1/2, 0), (1/2, 1) and (1/2,2), with 
positive parity, leading to bound p{u))D* states, one of them identified as the 1)^(2460) ( I{J^) = 
^(1~) ). Therefore, the resonance 1)2(2460) was generated as a pD* quasibound state or molecule 
by the strong and attractive pD* interaction. As discussed in Ref. [23;], because of the large binding 
energy per p meson in spin 2, it is possible to obtain bound systems with several p mesons as building 
blocks. As mentioned in Ref. [37|, the pD* interaction is also very strong and can bind the system. 
Thus the main aim in the present work is, first, to study the three-body interaction of D* and two 
p mesons, for which we have two options, the clusters D* — f2{1270){pp) and p — (2460) (/9-D* ) , in 
order to see if there are some resonance structures in the scattering amplitudes. If this is the case, 
one can predict a not-yet-discovered -D3 resonance, and we could continue our study extending 
these ideas to include more p mesons as building blocks of the many-body system. Then we repeat 
the test in the four-body system and so forth. In analogy to the K*-umlti-p systems we find 
clear peaks in the scattering amplitudes for system with increasing number of p mesons. Thus we 
predict new resonances in our work which are basically quasibound states or molecules made of an 
increasing number of p{770) mesons and one D* meson which are not yet reported in the list of 
the PDG [H. 
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II. pp AND pD* TWO-BODY INTERACTIONS 



To evaluate the Faddeev equations under the Fixed Centre Approximation, we need to define 
the two-body cluster and then let the third particle collide with the cluster. Thus the starting point 
of our work is the two-body pp and pD* interactions, which were studied in Refs. 36|] and 3711 with 



the local hidden gauge formalism 33-35] and the unitary coupled channels method 0-13 



151-UJ. 



36( 1 and 37[ here to explain how to obtain the unitarized 



We briefly summarize the model of Refs. 
pp and pD* scattering amplitudes. 

To evaluate the pp and pD* scattering amplitudes with the coupled channels unitary approach, 
the Bethe-Salpeter equation in coupled channels is used: 



t=[l- VG]-^V, 



(1) 



where the kernel y is a matrix of the interaction potentials in each channel, which is calculated 
from the hidden gauge Lagrangian. G is a diagonal matrix of the loop function of every channel. 
The details can be seen in Refs. HQ. 

To construct the three-body system we start from the clusters /2(1270) ( I{J^^) = 0(2"'"'") ) 
and L»^(2460) ( /(J^) = 1(1") ) and add to them a D* or a p respectively. The new particles are 
introduced with their spins aligned with that of the cluster such that the total spin adds one unity. 
Thus we only need to take into account the potential of spin 5 = 2 for pp and pD* interactions. 
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(2) 
(3) 
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where g = Mv/^/t^ with My the vector meson mass and /tt the pion decay constant. In these 
equations ki, i = 1,2,3,4 are the initial, (1,2), and final (3,4) momenta of the particles. The 
quantity k = rrip/m^, appears because in some transitions one is exchanging a heavy vector 



instead of a light one. Note that in isospin / 
2 is ojD*. 

As mentioned in Refs 



1/2 there are two coupled channels, 1 is pD* and 



36( 1 and 371], we also should take into account the contribution of the 



box diagram with two pseudoscalar mesons in the intermediate state. We only add the imaginary 
part of the box diagram contribution to the potential V, which is not accounted for by the coupled 
channels 3^, 40], and neglect the real part which is very small. We finally note that we take into 
account the p mass distribution by replacing the G function by its convoluted form by the mass 
distribution of the p mesons in the loop. 

Finally, considering the box diagram contribution to the potential V and the convolution of 
the p mass distribution in the loop function G, we show the evaluated results of pp and pD* 
scattering amplitudes in Fig. (H which are consistent with Refs. [s^ and 37]. The structure of the 
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FIG. 1: Modulus squared of the scattering amplitudes. Left: \ti=°\'^, /a (1270); Right: ^2(2460) 




FIG. 2: Diagrammatic representation of the FCA to Faddeev equations. 



resonances /2(1270) and 1)2(2460) are clear in the peak of the modulus squared of the amplitudes. 
The nonresonant amplitudes t^/p~^'^~'^^ and t 



{I=3/2,S=2) 
pD' 



are not shown here. 



III. MULTI-BODY INTERACTION FORMALISM 



The Faddeev equations under the Fixed Centre Approximation are an effective tool to deal 
with multi-hadron interaction 2314291]. They are particularly suited to study systems in which a 
pair of particles cluster together and the cluster is not much modified by the third particle. The 
Fixed Centre Approximation to Faddeev equations assumes a pair of particles (1 and 2) forming 
a cluster. Then particle 3 interacts with the components of the cluster, undergoing all possible 
multiple scattering with those components. This is depicted in Fig. [5J With this basic idea of the 
Fixed Centre Approximation, we could write the Faddeev equations easily. For this one defines 
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two partition functions Ti , T2 which sum ah diagrams of the series of Fig. [2] which begin with 
the interaction of particle 3 with particle 1 of the cluster (Ti), or with the particle 2 {T2). The 
equations then read 



Ti — ti + tiGoT2, 

T2 = t2 + t2GQTi, 

T = Ti + T2, 



(8) 
(9) 
(10) 



where T is the total three-body scattering amplitude that we are looking for. The amplitudes ti 
and t2 represent the unitary scattering amplitudes with coupled channels for the interactions of 
particle 3 with particle 1 and 2, respectively. Besides, Go is the propagator of particle 3 between 
the components of the two-body system. 

For the unitary amplitudes corresponding to single-scattering contribution, one must take into 
account the isospin structure of the cluster and write the ti and t2 amplitudes in terms of the 
isospin amplitudes of the (3,1) and (3,2) systems. Details will be discussed in next section. Besides, 
because of the normalization of Mandl and Shaw 4]| which has different weight factors for the 
particle fields, we must take into account how these factors appear in the single scattering and 
double scattering and in the total amplitude. This is easy and is done in detail in 24, [2^. We 
show below the details for the present case of a meson cluster (also particle 1 and 2) and a meson 
as scattering particle (the third particle). In this case, following the field normalization of 41[ we 
find for the S matrix of single scattering. 



SI 



(1) 




(11) 



(12) 



X {2TTr5{k + kR 



where, k, k' {kji, k'^) refer to the momentum of initial, final scattering particle {R for the cluster), 
V is the volume of the box where the states are normalized to unity and the subscripts 1, 2 refer 
to scattering with particle 1 or 2 of the cluster. 

The double scattering diagram, Fig. [2] (b), is given by 

1111111 



5(2) 



i{2TT)^S{k + kR 



k' - k'pi) 



R)y2 

1 



(13) 



,02 



ml + ie 



tlt2, 



t -1' 

where FR{q) is the cluster form factor that we shall discuss below. Similarly the full S matrix for 
scattering of particle 3 with the cluster will be given by 

1 



S = -iT {2Tr)^6{k + kR - k' - k'^) 
1111 
\/2a^ \/2^ ^2^' 



V2 



In view of the different normalization of these terms by comparing Eqs. (llip . (fT2 
(jl4p . we can introduce suitable factors in the elementary amplitudes. 



2mR 
2mi 



t2 



2mR 
2m2 



t2, 



(14) 



13]) and 



(15) 
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where we have taken the approximations, suitable for bound states, = ' ^'^'^ ^ 
diagrams by means of 

The function Gq in Eq. (fT6|) is given by 

J (zvr)'^ q'J^ — — mg + z e 

where Fji{q) is the form factor of the cluster of particles 1 and 2. We must use the form factor of 
the cluster consistently with the theory used to generate the cluster as a dynamically generated 
resonance. This requires to extend into wave functions the formalism of the chiral unitary approach 
developed for scattering amplitudes. This work has been done in 43-44| for s-wave bound states, 



s-wave resonant states and states with arbitrary angular momentum respectively, here we are only 
need the expressions for s-wave bound states, and then the expression for the form factors is given 



in section 4 of [43(], which we use in the present work and reproduce below 

F (q) - ^ [ d^p ^ ^ 

J\p\<A'\p-q]<A' 2Ei{p)2E2 



b1<AMp-,1<A' '^2E,{p)2E2{p}MR-E,{p)-E2{p} 
1 1 1 



2Ei{p- q) 2E2{p- q)MR- Ei{p- q) - E2{p - q) ' 

^ " Xi<A' '^'^(2i?i(p-) 2E2{p) Mr - E,{p) - E2{p)) ' ^^^^ 

where Ei and E2 are the energies of the particles 1, 2 and Mr the mass of the cluster. The 
parameter A' is a cut off that regularizes the integral of Eqs. (fT8|) and (fT9|) . This cut off is the 
same one needed in the regularization of the loop function of the two particle propagators in the 
study of the interaction of the two particles of the cluster [43|] . We take in the present work 
A' = 875 MeV, the same as used to generate the bound states |24l. [2^ of /2(1270) in the two-body 
interaction, and A' = 1200 MeV for getting the D2(2460) cluster. Thus we do not introduce any 
free parameters in the present procedure. 

In addition, q^ , the energy carried by particle 3 in the rest frame of the three particle system, 
is given by 

q'{s) = ^ - (20) 

Note also that the arguments of the amplitudes Tj(s) and ti{si) are different, where s is the 
total invariant mass of the three-body system, and Si are the invariant masses in the two-body 



systems. The value of Si is given by [24 1 



(M| + mf - mj){s -ml-Ml" 
2Ml 



mi + mi+'-^ ^,(.,J = 1,2, z^j) (21) 



where m/, (/ = 1,2,3) are the masses of the corresponding particles in the three-body system and 
Mr the mass of two body resonance or bound state (cluster). 
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IV. RESULTS 

The D*-multi-p interactions that we investigate in the present work are hsted in TablelU and are 
explained as follows. For the three-body interaction, we have two options: particle 3 = D* , cluster 
or resonance R = f2 (particle 1 = p, 2 = p) and 3 = p, R = D2 {I = p, 2 = D*). For four-body, 
we also have two cases: 3 = f2, R = D2 {1 = p, 2 = D*) and 3 = D2-, i? = /2 (1 = p, 2 = p). For 
five-body, 3 = D\ R = {I = /z, 2 = /a) and 3 = p, i? = Z)| (1 = /a, 2 = Dl). For six-body, 
2, = Dl, R = {I = /a, 2 = /s) and 3 = f2, R = 0% {I = /a, 2 = Dl). We describe all these 
cases in detail below. 



TABLE I: The cases considered in the Z?*-multi-p interactions. 



particles: 


3 


R (1,2) 


amplitudes 


Two- body 


P 


D* 


tpD' 




P 


P 


tpp 


Three-body 


D* 


/2 (PP) 


TD'-f2 




P 


D*2 (pD*) 


Tp-D; 


Four-body 


D* 








f2 


D*2 {pD*) 




Five-body 


D* 


h {f2h) 






P 


Dl {hD*2) 




Six-body 




h (/2/2) 






/2 


Dl {f2Dl) 





A. Three-body interaction 

For three-body interaction, we have two options of structure: D* — f2{pp) and p — D2{pD*), 
which means 3 = D* , R = f2 {1 = p, 2 = p) and 3 = p, R = {1 = p, 2 = D*). Thus, to 
evaluate these scattering amplitudes, we need as input the ti and t2 amplitudes of the (3,1) and 
(3,2) systems, ti = ^2 = ipD* for D* — f2{pp) and ti = tpp, t2 = tpD* for p — D2{pD*). We should 
calculate the two-body pp and pD* amplitudes. 

As mentioned before, the isospin structure of the cluster should be considered for the ti and t2 
amplitudes. For the case of D* — f2ipp), the cluster of /a has isospin / = 0. Therefore the two p 
mesons are in an / = state, and we have 

\pp >(0'0)= > > -1(0,0) > ), (22) 

where |(1,— 1) > denote |(/^,/^) > which shows the Iz components of particles 1 and 2, and 
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\pp means \pp >^^'^'\ The third particle is a D* meson taken |/|) >= |i >. Then we obtain 
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where the notation of the states followed in the terms is |(|,|),— 1 >= \ , I^^), l"^ > for tsi, and 
> for t32. Then we find 



(24) 



But for the case of p — D2{pD*), the situation is different. Because the isospins of p and D2 



are Ip = 1 and = ^, the total isospin of the three-body system are Itotai = IppD* 
Itotai = IppD* = |, and then we have 



or 



\pD* >(i2)= \ppD* 



\PD*2 > 



\ppD* > 



1(1, 
1(1, 
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Jl(o,^)>, 



Ix /2,. 1. 
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(25) 



where we have taken = ^ for convenience. Therefore the \pD* > states inside the D2 for the 
Iz = and Iz = +\ are given by 
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For the two possibilities, combining Eqs. (|25p and (|26p and performing a similar derivation of Eq. 
(I23l). we obtain 
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We show our results in Fig. [3l In Fig. [3] (left) we show the modulus squared of the amplitudes for 
IT^tIjjP ^^"^ l^p-r)*^P' ^^'^ ^^'^ there are clear peaks around the energy 2800 — 2850 MeV 
which is about 400 MeV lower than the D* — f2 threshold. The bindings are large because they scale 
with the mass of the mesons and we have now a D* interacting with two p mesons. The strength of 
the peak of {T^^^^J"^ is two times bigger than for \T^~^^\'^, and we see that the D* — f2 component 
is a bit more bound than the p — D2 one. We expect that a real state would be an admixture 
of both with a binding in between that of the individual components. In Fig. [3] (right) we show 
|rf~|l'l^P, and there is a clear resonant structure about 3120 MeV, the strength of which is 30 

times smaller than that of \T^~]it\'^ in the left figure and less bound. We are concerned with 
the lowest lying states and hence we concentrate on the predicted new state with a structure 
formed by a mixture of D* — f2 and p — D2, with a mass about 2800 — 2850 MeV and a width 
about 60 - 100 MeV. 



B. Four-body interaction 

There are also two possibilities in the four-body interaction as we have shown in Table |IJ particle 
3 = /2, i? = (1 = p, 2 = D*) or particle 3 = D^, R = f2 {1 = p, 2 = p). Because //^ = and 
lo* = ^, the total isospin of the four-body system is only Itotai = f • the first case, /2 collides 
with the L>2, the amplitudes ti = tf^p = Tp^f^ has been evaluated in Ref. [23] and is reproduced 
in our work, and ^2 = if2D* = ^D*-/2 which has been evaluated in the former subsection IIV Al For 
the second case, D2 collides with the /2, and the amplitudes ti = t2 = t^^p = Tp_£,* have been 
evaluated in the former subsection IIV Al We must now consider that the three-body amplitude 
Tp-D^ is also combined with different isospins as mentioned in subsection IIV Al This situation is 
similar to the case when the D* collides with the /2, because the isospins of both the D2 and D* 
are I = \, thus from Eq. we have 

ti = Tpm = \ {2tIT"' + TiT'") , t2 = t,. (28) 

The results are shown in Fig. HI The left of Fig. [5] is |T^7^j^p. We find that there is a clear 
peak at an energy of 3200 MeV, the width of which is about 200 MeV. The right of Fig. U] shows 
l^/2--D*l'^ and there is a resonant peak around the energy 3075 MeV with a large width of nearly 
400 MeV. The strength of the peak of \Ti~_}0, p is about two times bigger than the one of \T^^^? p 
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FIG. 4: Modulus squared of the TD*-f2 (left) and Tf^-D* (right) scattering ampHtudes. 



and the energy of the peak is more bound too. But from the former results, subsection IIV Al we 
found that iT^T^j^p has more strength and is more bound than iTj^"""^?]^. We have investigated 

that this is because of the contribution of |T'f~n'*^P, even though the strength of |T^~il{^P is much 

smaller and less bound than the one of \Ti=l(h'' fr om the former results. When we removed the 

contribution of in Eq. (|28p . the strength of the peak was enhanced by a factor five and 

was more bound. Therefore, within the uncertainty of the theory, we find a new D^^ resonance, of 
a mass about 3075 - 3200 MeV and a width about 200 - 400 MeV. 



C. Five-body interaction 



For the five-body interaction, we also have two options for the cluster, one of which is the 
particle /4 studied in Ref. [2^] and the other one the resonance obtained in the four-body 
interaction, subsection IIV B[ Thus letting the third particle {D* or p) collide with them, we have 
^ = D*, R = fi {1 = /2, 2 = /s) or 3 = p, i? = (1 = /2, 2 = D*). Because the isospin If^ = 
and Id* = i, the total isospin of the five-body system is only Itotai = | iii the D* — structure, 
but Itotai = I or hotai = | the p — D\ structure. Thus the situation is similar to the three-body 
interaction discussed before, D* (or p) collide with /2 (or 02)- Therefore in the first case, the 
D* collides with the /4, and the amplitudes ti = t2 = i_D*/2 = T^r}j:^^ have been evaluated in 
subsection IIV A[ For the second case, the p collides with the I?4, which is similar to p — D2 in the 
three-body interaction, thus, after doing a similar derivation as in Eq. (|23|) . we have 

^(/=l/2) . _ ^(/=1) _ ^/=l/2. 

t(^=3/2) . f _ 7.(/=l) f _f _ rpI=3/2 ^ ^ 

^p-Dl ■ *1 - V/2 - ^31 ' *2 - tpD* - J32 



where the T^l is the same as Tp-f^ in the subsection IIV Bl reproducing the results of Ref. 23l |. 



J ]^ / 2 I 3/2 

and T and T ~J, have also been evaluated in subsection IIV Al 

p ^2 P ^2 

In Fig. Owe show our results. The left of Fig. [5] is |T^7lj^P and we observe a resonant peak 
around the energy 3375 MeV with a width of less than 200 MeV. The right of Fig. [5] is \TpZ^*\'^ 
and \TpZo*\'^ ■ We find that there is a resonant structure in |T^Z£)*^P at the energy 3360 MeV, the 
width of which is about 400 MeV, and the position is very close to the one of iT^Tly^p. But the 
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strength of IT J* I is one order smaller than the one iTnT 7P. For it: there is no resonant 

structure. Therefore, within uncertainties, we also find a new Dt^ resonance, with a mass about 
3360 - 3375 MeV and a width about 200 - 400 MeV. 



D. Six-body interaction 

Similarly to the five-body interaction, we also have two options of the cluster for the six-body 
interaction, the particle U studied in Ref. and the resonance D\ obtained in subsection IIV B[ 
Now letting a resonance {D2 or /2) be the third particle and collide with them, we have 3 = D2, 
i? = /4 (1 = /2, 2 = /a) or 3 = /2, i2 = Dl (1 = /a, 2 = D*). Because If^ = If, = and 
Id* = Idi = the total isospin of the six-body system is only Itotai = ^- Thus, in the first 
case, the D2 collides with the f^, the amplitudes ii = ^2 = iz)*/2 = TD*~-f2^ have been evaluated in 
subsection llVBI For the second case, the /a collides with the the amplitudes ti = tf^f^ = Tf^_f^ 



reproduce the results from Ref. [23|], and ^2 = tf-iD* = Tf2-D* has been calculated in subsection 
lYBl 

Our results are shown in Fig. El The left of Fig. [6] is \Tl^]!f P where we see a peak around the 

2 ■'4 

energy 3775 MeV with a large width of nearly 400 MeV. The right of Fig. [6] is iTf^^/,* P, and there 
we find that there is not a clear peak at the energy 3550 MeV. It looks like the resonant structure 
of /2 — D\ is not as stable as the ~ fi o'^^. From these results, we could predict a new 
resonance with more uncertainty, with a mass of about 3775 MeV and a width about 400 MeV. 



V. CONCLUSIONS 



In the present work, we show the results of our investigation of the D*-mu\t\-p systems. Our 
idea is based on the fact that the two-body interactions of pp and pD* in spin 5 = 2 are so strong 
as to bind the particles forming the resonances/2(1270) [32] and 1)2(2460) ^3] respectively. So 
we could study the many-body D*-um\ii-p systems in an iterative way looking at the structure of 
the amplitudes and observing clear peaks that become wider as the number of p mesons increase. 
The work proceeded analogously to the study of multi-/) system in [ill, where the p3(1690)(3 ), 
/4 (2050) (4++), p5(2350)(5 ), and /e (2510) (6"*"+) were described as basically molecules of multi- 
p(770) states, and similarly to the work of Q where the ir^(1430), K^(1780), i^|(2045), i^^(2380) 
and -fCg could be interpreted as molecules made of one i^*(892) meson and an increasing number 
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of p{770) mesons. The D*-multi-p states with spins ahgned combined to give some new charmed 
resonances, -D3, D^, and Dg, which are basically made of one D* meson and an increasing 
number of p(770) mesons and are not found in the list of PDG [38|. Their masses are predicted 
around 2800 - 2850 MeV, 3075 - 3200 MeV, 3360 - 3375 MeV and 3775 MeV respectively. And 
their widths are about 60 - 100 MeV, 200 - 400 MeV, 200 - 400 MeV and 400 MeV respectively. 
The analogy with the states already known in the strange and non-strange sector, together with 
the stronger interaction of the D* mesons, make our predictions solid within the uncertainties 
admitted. We are, thus, reasonably confident that such states can be found in the future in coming 
facilities like FAIR and others. 
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